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We simulate direct photon production in evolution of viscous QGP medium. Photons from Comp- 
ton and annihilation processes are considered. Viscous effect on photon production is very strong and 
reliable simulation is possible only in a limited px range. For minimally viscous fluid r7/s=0.08), di- 
rect photons can be reliably computed only up to px < 1.3 GeV. With reduced viscosity (rj/s=0.04), 
the limit increases to px <2GeV. 
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I. INTRODUCTION 

Recent experiments at Relativistic Heavy Ion Col- 
lider (RHIC) have produced convincing evidences that in 
t/s nn =20Q GeV Au+Au collisions, a collective medium 
is produced [J-Q • Strong elliptic flow in non-central 
collisions is key evidence to this understanding. Hy- 
drodynamical analysis of experimental charged particles 
data, also suggests that in Au+Au collisions, a collec- 
tive medium, with viscosity to entropy ratio close to the 
AdS/CFT lower bound of viscosity, r]/s > 1/Att is pro- 
duced 0-0| • Whether the collective medium is a decon- 
fined Quark-Gluon Plasma (QGP) or not, is still a ques- 
tion of debate. Hopefully, the issue will be settled in 
Pb+Pb collisions at the Large Hadron Collider (LHC). 

Direct photons probe the early medium produced in a 
collision better than the charged hadrons. Hadrons, be- 
ing strongly interacting, are emitted from the surface of 
the thermalised matter and carry information about the 
freeze-out surface only. They are unaware of the condi- 
tion of the interior of the matter and can provide infor- 
mation about the deep interior only in an indirect way. In 
a hydrodynamic model, one fixes the initial conditions of 
the fluid such that the " experimental" freeze-out surface 
is correctly reproduced. In contrast to hadrons, photons, 
being weakly interacting, are emitted from whole volume 
of the matter. Throughout the evolution of the matter, 
photons are emitted. Conditions of the produced mat- 
ter, at its deep interior, are better probed by the pho- 
tons. Depending on the transverse momentum, direct 
photons can probe different aspects of heavy ion colli- 
sions. A thermalised medium of quarks and gluons, or of 
hadrons, can produce significant number of thermal pho- 
tons. They are low px photons (px < 3 GeV/c). Low 
Pt photons can test whether or not, QGP is produced in 
Au+Au collisions. Hard photons (px > 6 GeV/c) are of 
pQCD origin and test the pQCD models. Fast partons 
from 'jet' can interact with thermal partons of QGP and 
produce photons. At intermediate px range, (3 < px < 
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6 GeV/c), interaction jets with QGP could be an impor- 
tant source of direct photons [H [t| . 

In ideal hydrodynamic models, photon production in 
Au+Au collisions at RHIC energy has been studied ex- 
tensively [13, El- However, it is now realized that the 
strongly interacting medium, produced in Au+Au col- 
lisions, must be treated as a viscous medium. Gravity 
dual theories suggest that specific viscosity, i.e. viscosity 
to entropy ratio of any matter has a lower bound, the 
so called KSS bound rj/s = 1/4tt fl3,ll3|. Even though, 
photons are important probe of QGP matter, viscous ef- 
fects on photon production are not studied much. Only 
recently, Dusling [l4[ studied viscous effects on photon 
production from a QGP medium. However, the model 
appears to have some inconsistency. In viscous evo- 
lution, photon production is affected due to (i) modi- 
fied fluid evolution and (ii) non-equilibrium correction 
to equilibrium distribution function. In [14}, while the 
non-equilibrium correction to the distribution function 
was included, modification of the fluid evolution, due to 
viscosity was neglected. In the present paper, the incon- 
sistency is removed. 



II. HYDRODYNICAL EQUATIONS, EQUATION 
OF STATE AND INITIAL CONDITIONS 

In a hydrodynamical model, the invariant distribution 
of direct photons is obtained by convoluting the photon 
production rate with space-time evolution of the fluid. 
We assume that in ^^^=200 GeV, Au+Au collisions 
at RHIC, a baryon free QGP fluid is formed. Space-time 
evolution of the fluid is obtained by solving 2nd order 
Israel-Stewart's theory, 

d^ v = 0, (1) 

Dtt^ = - — (n" 1 ' - 2r]V <fi u v> ) 

- [yV A +u u it yX ]Du x . (2) 

Eq[T] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v — pg^ v + , 
e, p and u being the energy density, pressure and fluid 
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velocity respectively. 7r M " is the shear stress tensor 
(we are neglecting bulk viscosity). EqJ5] is the relax- 
ation equation for the shear stress tensor tt^ . In Eq|2j 
D = u^dfj, is the convective time derivative, V <m m 1/> = 
±(V^ + Vu") - \{d.u){g^ ~ u^u") is a symmetric 
traceless tensor. 77 is the shear viscosity and r v is the 
relaxation time. It may be mentioned that in a confor- 
mally symmetric fluid relaxation equation can contain 
additional terms [l5j . Assuming boost- invariance, EqsJT] 
and [5] are solved in (r = \/t 2 — z 2 , x, y, rj s = ^lnf^f) 
coordinates, with the code " 'AZHYDRO-KOLKATA" ', 
developed at the Cyclotron Centre, Kolkata. Details of 
the code can be found in [f| 0, EH ■ 

Eqs ll|2l are closed with an equation of state p = 
p{e). Lattice simulations pL7H20| indicate that the 
confinement-dcconfincment transition is a cross over, 
rather than a 1st or 2nd order phase transition. There is 
no critical temperature for a cross-over transition. How- 
ever, one can define a pseudo critical temperature, the 
inflection point on the Polyakov loop. In Wuppertal- 
Budapest simulation [l9l.[20j. pseudo critical temperature 
T c «170 MeV. In the present simulation, for the QGP 
phase, we use an EOS based on Wuppcrtal-Budapest sim- 
ulation. 

Solution of partial differential equations fEqs lll2j) re- 
quires initial conditions, e.g. transverse profile of the 
energy density (e(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (tt^(x, y)) at the initial time Tj. 
One also need to specify the viscosity (f?) and the relax- 
ation time (r w ). A freeze-out temperature is also needed. 
In the following, we will consider viscous effects on pho- 
tons from the QGP phase only. The hydrodynamical 



equations are then solved till the freeze-out temperature 
Tp = T c =170 MeV. At the initial time tj, initial energy 
density is assumed to be distributed as @ 



e(h,x,y) = Ei{(l-f hard N part (b,x,y)+f hard N coU (b,x,y)], 

(3) 

where b is the impact parameter of the collision. N part 
and N co u are the transverse profile of the average partic- 
ipant and collision number respectively, calculated in a 
Glauber model. fhard=0-13 is the hard scattering frac- 
tion [2l|. Si is the central energy density of the fluid 
in impact parameter b = collision. As it will be dis- 
cussed later, we have simulated Au+Au collisions for 
a range of initial (central) energy density and initial 
time. We also assume that initial fluid velocity is zero, 
v x (x,y) = v y (x,y) — 0. The shear stress tensor was ini- 
tialized with boost- invariant value, tt xx = n vv = 2i]/3ri, 
ir xy =0. For the relaxation time, we use the Boltzmann 
estimate r n — 3r//Ap. We also assume that shear viscos- 
ity to entropy ratio is a constant throughout the evolu- 
tion. In the following, we simulate Au+Au collisions for 
?7/s=0-0.12. 

III. PHOTON RATES 

As mentioned earlier, in viscous evolution photon rates 
are modified. The photon rate equations involve distri- 
bution functions of quarks and gluons. For example, in 
1 + 2 — > 3 + 4 processes (e.g. Compton and annihilation 
processes), the general form of photon production rate 
is, 



dR 

Ei—=N 



d?p l d 3 p 2 
(2tt) 3 2E 1 (2it) 3 2E 2 



h(El)f 2 (E 2 )(2ir)6 4 ( Pl +P2-P3- P4)|M| 2 [1 ± f 3 (E 3 )] 



d 3 P 3 



(2ir) 3 2E 3 



(4) 



where M is the matrix element for the reaction, /'s in 
EqH]are distribution function of quarks and gluons. Un- 
like an ideal fluid, in viscous evolution, each distribution 
functions / is modified due to non-equilibrium correction, 



7ieq 



f eq (l + Sf neq ), 5f neq « 1. 



(5) 



The non-equilibrium correction 6f neq depend on dissi- 
pative forces as well as on particle momenta. For shear 
viscosity, it can be obtained as, 



Sfneq = Cp^TT^ 



1 



(6) 



2{e+p)T 2 

It is obvious that non-equilibrium correction to photon 
rates is non-trivial. For Compton and annihilation pro- 
cesses, in the leading log approximation, the rate equa- 
tion can be simplified 14], 



scr(s) 
Eo 



(7) 



Non-equilibrium correction to distribution functions 
inside the phase space integral leads to corrections in the 
logarithmic term and can be neglected. For equilibrium 
distribution function, as demonstrated by Kapusta and 
Lichard [22[ , the phase integration can be explicitly eval- 
uated and one obtain for both Compton and annihilation 
processes 



,dR 
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E^r- fa ^tf nea (E)T z In 

d 3 p 9 2ir 2 J qy ' 
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The invariant photon distribution then has two parts, 



equilibrium part (E d ^" q ) and a non-equilibrium part 



3 



10 1 
10° 

f 10- 2 
i— 

10- 4 

>> 

z 10- 5 

10 6 

10- 7 

1 2 3 4 5 

PT (GeV) 

FIG. 1: (color online) Transverse momentum spectra for 
photons, from evolution of ideal and minimally viscous 
(?7/s=0.08) QGP. For viscous fluid, spectra obtained with and 
without the non-equilibrium correction to distribution func- 
tion, are shown separately. 



(E dI ^3° q )■ Since non-equilibrium correction to distribu- 
tion function is assumed to be small, it is essential that, 
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IV. 



EFFECT OF VISCOSITY ON PHOTON 
SPECTRA AND ELLIPTIC FLOW 



Gravity dual theories [12| indicate that viscosity to 
entropy ratio of strongly interacting matter is bounded 
from the lower side, r\j s > \/An. We first consider ef- 
fect of minimal viscosity on photon production. We as- 
sume that minimally viscous (rj/s—0. 08) QGP fluid is 
thermalised in the time scale Ti=0.6 fm to initial central 
temperature Tn=350 MeV. A large number of charged 
particle's data e.g. identified particles spectra, elliptic 
flow etc. data are explained in hydrodynamical model 
with similar initial time and temperature scale [f|. In 
FigQ] simulated photon spectra in 20-40% Au+Au col- 
lisions, from evolution of ideal and minimally viscous 
QGP fluid, are shown. The spectra with or without 
the non-equilibrium correction are shown separately. If 
non-equilibrium correction to distribution function is ne- 
glected, pt spectra of photons is increased by a factor 
of ~ 1.2-1.5. The increase is largely pt independent. 
In contrast, when non-equilibrium correction is included, 
photon production is increased more at large px than at 
low pt, It is also expected, non-equilibrium correction 
increases with px (see Eq[S]). The arrow in Fig Q] indicate 
the approximate px when non-equilibrium contribution 
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FIG. 2: (color online) same as in Fig[T]but for elliptic flow. 



(SN neq ) to photon spectra equals the equilibrium contri- 
bution N eq ). For minimally viscous fluid, the equality 
occur at pr ~1.7 GeV. As noted earlier, viscous hydro- 
dynamics is applicable only when SN neq /N eq << 1. Ev- 
idently, for ?7/s=0.08, in a hydrodynamic model, photon 
production from Compton and annihilation processes can 
not be reliably computed beyond pr=l-7 GeV. Indeed, 
if we assume that viscous hydrodynamics remain reli- 
able only until, S ^" q = 0.5, the photons from Compton 
and annihilation processes can be computed only up to 
p T -1.3 GeV. 

In Fig 13J we have shown the simulation results for el- 
liptic flow for photons. At very low px elliptic flow is 
negative. At low px < 0.5GeV present model is not 
reliable. The logarithmic factor in the photon rate equa- 
tion is negative when E < 0.27g 2 T. For fluid (central) 
temperature Tn=350 MeV, reliable computation is pos- 
sible only beyond px —0.5 GeV. In ideal fluid evolution, 
at px >0.5 GeV, elliptic flow increases with px, till a 
maxima is reached at px ~2 GeV. At even higher px, it 
decreases again. Note that even the highest V2(px) is not 
large, less that ~ 2%. Compare this value to v 2 (px) ~ 
20% for charged particles. In viscous evolution, elliptic 
flow is reduced further. Incidentally, very small photon 
elliptic flow is consistent with experiments. In experi- 
ments also, photons do not flow |27], [28j]. In Fig 121 one 
observes that if non-equilibrium correction is neglected, 
elliptic flow in viscous evolution is more than flow in ideal 
fluid evolution. However, with non-equilibrium correc- 
tion included, flow is reduced in viscous evolution. It 
shows that it is important to have a consistent model, 
otherwise, one can conclude wrongly. The arrow in Figj2] 



eg 



N, 



indicate that approximate px when SN n 

Photon spectra from Compton and annihilation pro- 
cesses, for four values of viscosity to entropy ratio, r]/s= 
(ideal fluid), 0.04, 0.08 (AdS/CFT lower bound) and 0.12 
are shown in FigOl Initial time and temperatures are 
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FIG. 3: (color online) Photon spectra from evolution of QGP 
fluid, with viscosity to entropy ratio, r//s=0, 0.04, 0.08 and 
0.12 are shown. The initial time is n=0.6 fm and initial cen- 
tral temperature is Ti=350 MeV. 



Ti=0.6 fm and Ti=350 MeV. As expected, high px yield 
is increased in more viscous fluid evolution. The arrows 
in Fig (3] indicate the approximate pr when SN neq = N eq . 
Photon production can not be computed reliably beyond 
p T =l-3, 1.7 and 2.4 GeV, for QGP viscosity, r?/s=0.12, 
0.08 and 0.04 respectively. Very limited pr range over 
which viscous hydrodynamics remain applicable for pho- 
ton production raises certain issues, which will be dis- 
cussed later. We just mention that viscous effect on pho- 
ton production is much stronger than in charged par- 
ticles. Indeed, for charged particles, viscous hydrody- 
namics remain applicable over a much wider pt range 
[S 13 ■ The reason is understood. Photons are emitted 
through out the evolution. At early time, shear stress 
tensors have finite values and non-equilibrium correction 
to photon production is large. In contrast, charged parti- 
cles are emitted from the freeze-out surface. Shear stress 
tensors evolves very fast. Compared to early times, stress 
tensors at freeze-out are much smaller in magnitude. Ac- 
cordingly, non-equilibrium correction is small. We have 
not shown simulation results for elliptic flow as a func- 
tion of viscosity. As shown earlier, V2(j>t) for photons 
is very small in ideal fluid evolution. It further reduces 
with r)/s. 

Since photons are emitted throughout the evolution, 
including the earliest phase of QGP, one hope to ex- 
tract QGP formation time from photon measurements. 
Indeed, in [25|, [26| it was suggested that photon ellip- 
tic flow can be used to constrain QGP formation time. 
However, from the experiments [U Ha], it appears that 
photons do not seem to experience flow. One can not 
possibly extract QGP formation time from elliptic flow 
measurements. In 14], Dusling suggested that inverse 
slope of the photon spectra, if measured accurately, can 
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FIG. 4: (color online) Initial time dependence of photon spec- 
tra from evolution of QGP fluid with viscosity 7//s=0.08. The 
initial entropy density times the initial time is fixed at SiTi=60 
fm~ 2 . 



be used to put stringent bound on QGP viscosity and 
initial (formation) time. 

In Fig 01 simulated photon spectra from evolution of 
minimally viscous (?7/s=0.08) QGP fluid, for four values 
of initial time Tj=0.2, 0.6, 1.0 and 1.4 fm are shown. The 
initial temperature is obtained from the condition that 
initial entropy density times the initial time is a constant, 
SiTi — 60/m" 3 . As the initial time is reduced, photon 
spectra are hardened. Photon yield is also increased. It is 
easily understood, for fixed SjTj, initial time and temper- 
ature are inversely related. As the initial time is reduced, 
high px photon production is increased due to increased 
fluid temperature. The arrows in FigJH indicate the pr 
when non-equilibrium correction to spectra is equal to 
the equilibrium contribution. 

To obtain the inverse slope parameter T e tf, we have fit- 
ted the spectra in the pr range 1.5 < pr < 2.5GeV with 
an exponential (dN/cPpx oc e~ PT / T ' SI ). Ideal hydrody- 
namic simulations for photon spectra suggest that in this 
Pt range, QGP photons dominate the spectra [ToL [i"lj]. 
In Figl5l for fluid viscosity, rj/s—0, 0.04 and 0.08, inverse 
slope parameter is shown as a function of initial time. 
Inverse slope parameter T e ff decreases with increasing 
initial time, as well as with decreasing viscosity. In FigJSJ 
the shaded region indicates the experimental slope mea- 
sured in the PHENIX experiment. For ideal fluid, sim- 
ulated T e ff agree with experiment for r, > 1 fm. For 
viscous fluid, in the time scale, 0.2 — 1.0 fm, inverse slope 
of the simulated spectra are higher than that observed in 
experiment. However, we have neglected photons from 
the hadronic phase. Hadronic photons will dominantly 
contribute at low pt, reducing the inverse slope param- 
eter. In other word, in more realistic simulation, T e ff 
could be smaller than obtained presently and even for 
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FIG. 5: (color online) Effective temperature from photon 
spectra for initial time Ti=0.2, 0.6, 1.0 and 1.4 fm. The ini- 
tial entropy density times the initial time is fixed at s;T; = 60 
fm~ 2 . Results are shown for viscosity to entropy ratio r//s=0, 
0.04 and 0.8 respectively. The two dashed lines indicate the 
experimental uncertainty in PHENIX measurements. 



viscous fluid, depending on the initial time, applicabil- 
ity range is less than pt=1-5 GeV. Applicability range is 
even less in more viscous fluid. Present simulations ne- 
glect Bremsstrahlung photons in the QGP phase. How- 
ever, their inclusion will not improve the situation. One 
wonders about the photons in the pr range 1.5 < Pt < 3 
GeV. As noted earlier, ideal hydrodynamics simulations 
indicate that in this px range, thermal photons domi- 
nate. However, present simulations indicate that if QGP 
fluid is minimally viscous, the photons in the px range 
1.5 < pt < 3 GeV can not be considered thermal and can 
not be described hydrodynamically. Do we understand 
them as pQCD photons? Or are they from a fluid with 
viscosity much less than the AdS/CFT lower bound? In- 
deed, in certain gravity dual theories, KSS bound can 
be violated by ~36% [29[. Limited pt range over which 
photons can be reliably computed in viscous dynamics 
also indicate that within a same model, one possibly can 
not explain, charged particle production and photon pro- 
duction. 



SUMMARY 



viscous fluid could be in agreement with experiment. 

To constrain initial time and viscosity from inverse 
slope of photon spectra will not be an easy task. As 
seen in FigjSJ inverse slope parameter from evolution of 
low viscosity fluid initialized at small Tj could be confused 
with T e ft from high viscosity fluid initialised at large Tj . 
As shown in Figj5l for rj/s = 0.04±0.04, depending on the 
initial time, change in inverse slope parameters is ^20-80 
MeV. Depending on the realistic situation, measurement 
of inverse slope parameter with accuracy 5T e f f rs 20-80 
MeV could possibly give a range of values for initial time 
and viscosity. 

Before we summarise, few comments are in order. We 
have shown that effect of viscosity is quite strong on 
photon production. Indeed, viscous hydrodynamics in 
applicable only in a limited px range. For minimally 



To summarize, we have studied effect of shear viscos- 
ity on Compton and annihilation photons. In viscous dy- 
namics, photon production is modified due to (i) changed 
space-time evolution of the fluid and (ii) non-equilibrium 
correction to the equilibrium distribution function. The 
non-equilibrium correction grows with viscosity as well 
with transverse momentum. Viscous effects on photon 
production are strong. Even for AdS/CFT lower bound 
of viscosity (rj/s=0.08), strong viscous correction render 
the hydrodynamics inapplicable beyond pr ~1.5 GeV. If 
QGP viscosity is larger than the Ads/CFT limit, then 
limitation will be even more. Photon production as 
a function of initial time, also suggest that if the in- 
verse slope parameter of the photon spectra, is measured 
within an accuracy of ±20 MeV, one can possibly limit 
the initial time and viscosity. 
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